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Background: COPD is a heterogeneous disease characterized by airflow obstruction and diag- 
nosed by lung function. CT imaging is emerging as an important, noninvasive tool in phenotyping 
COPD. However, the use of CT imaging in defining the disease heterogeneity above lung func- 
tion is not fully known. 

Methods: Seventy-five patients with COPD (58 men, 17 women) were studied with CT imaging and 
with measures of airway inflammation. Airway physiology and health status were also determined. 
Results: The presence of emphysema (EM), bronchiectasis (BE), and bronchial wall thickening 
(BWT) was found in 67%, 27%, and 27% of subjects, respectively. The presence of EM was associ- 
ated with lower lung function (mean difference % FEV^, -20%; 95% CI, -28 to -11; P< .001). 
There was no difference in airway inflammation, exacerbation frequency, or bacterial load in 
patients with EM alone or with BE and/or BWT ± EM. The diffusing capacity of the lung for car- 
bon monoxide/alveolar volume ratio was the most sensitive and specific parameter in identifying 
EM (area under the receiver operator characteristic curve, 0.87; 95% CI, 0.79-0.96). Physiologic 
cluster analysis identified three clusters, two of which were EM predominant and the third char- 
acterized by a heterogeneous combination of EM and BE. 

Conclusions: The application of CT imaging can be useful as a tool in the multidimensional 
approach to phenotyping patients with COPD. CHEST 2011; 140(3) .634-642 

Abbreviations: BE = bronchiectasis; BWT = bronchial wall thickening; Dlco = diffusing capacity of the lung for carbon 
monoxide; EM = emphysema; RV = residual volume; Va = alveolar volume 



COPD is a heterogeneous disease characterized by 
airflow obstruction that is not fully reversible and 
is associated with a progressive decline in lung func- 
tion. ^ Airway inflammation in COPD is usually neu- 
trophilic,^'^ but eosinophilic airway inflammation has 
been found at stable state^^ and during exacerbations.^^ 
Pathologic processes in COPD include destruction of 
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lung parenchyma and changes within the large and 
small airways. s This evident disease heterogeneity 
cannot be defined by the FEV^ alone, and alternative 
methods need to be sought to develop mechanistic, 
prognostic, and therapeutic applications in the man- 
agement of COPD.^ CT imaging has emerged as a 
noninvasive tool in the "phenotyping" of COPD, with 
measures investigating changes in the airway wall and 
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lumen and within the lung parenchyma/^-^^ as well 
as assessments of the degree of emphysema (EM) 
and burden of small airways disease. CT imaging 
quantitative assessments have shown, inter alia, asso- 
ciations with parameters of airway physiology,!^-!^ 
important COPD outcome measures, ^^-^^ and sys- 
temic inflammatory mediators. 

However, the role of CT imaging in defining COPD 
phenotypes above and beyond conventional charac- 
terization with full lung function (spirometry, static 
lung volumes, and gas transfer) remains to be fully 
determined. Furthermore, the validity of CT imag- 
ing as a biomarker of airways disease has yet to be 
established. 

We hypothesized that CT imaging provides addi- 
tional value to clinical and physiologic parameters 
in the multidimensional phenotyping of COPD. We 
tested our hypothesis by evaluating whether pheno- 
typing of COPD using factor and cluster analysis 
with lung function parameters (physiologic clustering) 
would generate the radiologic COPD disease groups 
seen in clinical practice (EM, bronchial wall thick- 
ening [BWT], and bronchiectasis [BE]). 

Materials and Methods 

Patients 

Patients with a physician diagnosis of COPD as per GOLD 
(Global Initiative for Chronic Obstructive Lung Disease) criteria^ 
were recruited consecutively from general respiratory clinics and 
through local advertising to enter the longitudinal Biomarkers in 
COPD Exacerbation study. Patients with a diagnosis of asthma, 
current active pulmonary TB, or any other clinically relevant 
lung disease were excluded. All patients had obstructive spirom- 
etry with a postbronchodilator FEV/FVC ratio of < 0.7, whereas 
severity was classified according to the GOLD criteria. ^ Patients 
with COPD who demonstrated bronchodilator reversibility were 
not excluded. All patients gave written informed consent, and 
the study was approved by the Leicestershire, Northamptonshire, 
and Rutland ethics committee (REG 07/H0406/157). 

Study Design 

Seventy-five patients enrolled in the Biomarkers in COPD 
Exacerbation study, who had previously undergone a CT scan as 
part of their clinical management, were studied. The median 
(interquartile range) interval time between patient characteriza- 
tion and CT scanning was 15 (25) months. 

Measurements 

All patients had complete demographic data recorded, includ- 
ing age, duration of symptoms, and full smoking and medical 
history. Full lung function, including reversibility testing with 
400 |jLg inhaled albuterol, was also performed according to the 
American Thoracic Society/European Respiratory Society con- 
sensus guidelines. 22-26 Health status and symptom scores were mea- 
sured using the St. George Respiratory Questionnaire (University 
of London; London, England),^^ a two-part questionnaire measur- 
ing recollection of symptoms and assessment of current activity 
and impact, in which scores are expressed as a percentage of 



overall impairment, with zero indicating best possible health sta- 
tus and 100 indicating worst possible health status; the Chronic 
Respiratory Disease Interviewer-Administered Questionnaire 
(McMaster University; Hamilton, Ontario, Canada), which mea- 
sures disease-specific quality of life with scores ranging from 1 to 7, 
with a higher score indicating a better quality of life for the 
domains of emotion, mastery, dyspnea, and fatigue; and the 
visual analog scale (VAS) for the domains of cough, breathless- 
ness, sputum production, and sputum purulence, which consists of 
a 100-mm line for each measured domain, with "no symptoms" at one 
end and "the worst symptoms ever" at the other. ^9 Body composi- 
tion was calculated using BMI.^^ Spontaneous or induced sputum 
was collected and analyzed for bacteria,^^ including colony-forming 
units.32 Quantitative real-time polymerase chain reaction was used 
to estimate the total bacterial load based on the abundance of 
16S ribosomal subunit encoding genes. In brief, bacterial DNA 
was extracted from homogenized sputum according to the 
manufacturer's instructions, using the QIAmp DNA Mini Kit 
assay (QIAGEN, Ltd; Hilden, Germany). Quantification of the 
total bacterial load was performed using the SYBR green assay 
(PE Applied Biosystems; Warrington, England). Sputum samples 
were then processed to produce cytospins for cell differential. ^^'^^ 
Venous blood was collected for assessment of peripheral blood 
differential cell counts and serum C-reactive protein. 

CT Imaging 

CT imaging was performed using a Sensation 16-slice scanner 
(Siemens Healthcare; Knoxville, Tennessee). Sequential scanning 
was performed at maximal inspiration from the apex to the dia- 
phragm at 10-mm increments, with 1-mm collimation, while patients 
were in the supine position. Image reconstruction was per- 
formed by using a high-spatial-frequency algorithm through a 
512 X 512 matrix with a small field of view. Scanning time ranged 
from 30 to 45 s with a 120-kV peak and an effective tube current 
of 140 mA. Images were reported at a window width and level of 
1,600 and —500 Hounsfield units, respectively. CT scan qualita- 
tive radiologic reporting for EM, BE, and BWT was performed in 
accordance with published guidelines. BE was present when 
one or more of the following criteria were fulfilled: (1) an internal 
diameter of the bronchus greater than that of the adjacent pul- 
monary artery, (2) a lack of tapering of the bronchial lumen toward 
the periphery, or (3) visualization of the bronchus within 10 mm 
of the pleural space. EM was present when focal areas of low 
attenuation without visible walls were detected, and included the 
presence of centrilobular, panlobular, and paraseptal EM. BWT 
was assessed subjectively as described previously.^^ Previous stud- 
ies at our institution have shown that the qualitative description 
method of EM, BWT, and BE has good interobserver agreement 
and image reconstruction; scanning and reporting was optimized 
for analysis of BWT, EM, and BE as previously described. 

Statistical Analysis 

Statistical analysis was performed using PRISM, version 4 
(GraphPad Software; San Diego, California) and SPSS, version 16 
(SPSS Inc; Chicago, Illinois). Parametric and nonparametric data 
are presented as mean (SEM) and median (interquartile range) 
unless stated otherwise. Log-transformed data are presented as 
geometric mean (95% CI). For comparison of unpaired or paired, 
parametric or nonparametric groups, the Student t test, paired 
t test, Mann-Whitney test, and Wilcoxon matched pairs test 
were used, respectively. For comparison of three groups or more 
for parametric and nonparametric variables, the one-way analysis 
of variance or Kruskal-Wallis test was used and the test was 
used for proportions. Logistic regression analysis was used^^ to 
assess the relationship of the dependant variables of (1) presence 
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of EM, (2) presence of BE, and (3) presence of BWTwith explan- 
atory (independent) variables, using the block entry method. Vari- 
ables entered into the logistic regression model, chosen for clinical 
relevance, were the following continuous variables: FEV^ % pre- 
dicted, diffusing capacity of the lung for carbon monoxide (Dlco)/ 
alveolar volume (Va), residual volume (RV) % predicted, total 
sputum neutrophil count, and percentage sputum eosinophils (both 
log transformed). Model inspection showed no multicollinearity, 
goodness of fit was performed using the Hosmer-Lemeshow 
test, and the NagelKerke was used to estimate the variance 
explained by the model. The Wald test was used to assess the 
significance of the individual dependant variables to the model. 
Unsupervised multivariate modeling using principal component 
factor reduction analysis (orthogonal varimax rotation method) 
was used to explore airway physiologic pattern expression in the 
patients. Hierarchic cluster analysis was then applied from the 
identified factors to determine physiologic clusters of COPD. 
Clinical characteristics for the physiologic clusters were tabulated. 
One-way analysis of variance, Kruskal Wallis, and the test were 
used to compare parametric, nonparametric, and proportions of 
clinical characteristics between physiologic cluster groups. A P value 
of < .05 was deemed to be statistically significant. 



Results 

CT data were available in 75 patients (58 men, 
17 women). The most common clinical indication for 
CT scanning was investigation or exclusion of malig- 
nancy (40%), followed by assessment for radiographic 
evidence of BE (30%) and suitability for lung volume 
reduction surgery (20%). The clinical characteristics 
of patients are shown in Table 1. 

EM, BE, and BWT Are Common Radiologic 
Observations in COPD, With Significant Overlap 

The CT scan description of EM, BE, and BWT 
was present in 67%, 27%, and 27%, respectively. 
The number of patients with EM, BE, and BWT, and 
the associated overlap, are shown in Figure 1. No evi- 
dence of EM, BE, or BWT on CT scanning was 
found in 13 patients in the COPD cohort. Subjects 
with CT scan evidence of disease had worse lung 
function, but there was no difference in demo- 
graphic distribution, health status, or symptom scores 
(Table 2). 

Conventional Lung Function Reliably Identifies 
Patients With Radiologic EM and COPD 

CT scan evidence of EM compared with no evi- 
dence of EM was associated with lower lung func- 
tion (39 vs 61; mean difference, -22%; 95% CI, -30 
to -13; F<.001), impaired Dlco (3.8 vs 5.6; mean 
difference, -1.8; 95% CI, -2.7 to -0.9; P< .001), and 
increased airway obstruction (FEV^/FVC) (0.43 vs 0.57; 
mean difference, -0.14; 95% CI, -0.20 to -0.08; 
F< .001). The receiver operator characteristic curve 
(area under the curve) for the DlcoA^A to correctly 



Table 1 — Clinical Characteristics of Patients With 
COPD With Corresponding CT Scan 



Characteristic 


Patients With COPD With 
CT scan Data (n = 75) 


ividie, i\o. \ /c ) 


DO [i 1 ) 


Age, y, mean (range) 




VjUllCllL SlllUKci, iMU. \ /O I 


24 (32) 


JUjA-hiiiUKer, iMU. \ /c / 






AR (^0 1 CC'^) 


Exacerbation history in previous y 


4.3 (0.3) 


JjiVli, Kg/ 111 


oa 1 (f) a\ 


TPS iisap-p No (%) 


68 (91) 




1494 (Q9) 






FFV % nrprlioi-prla 


4fi (9) 


!-< /2i\7/2i-rc;i r^i nl-\7' 
±ic vci oiUiiiLy, /V 




FEV/FVC,^ % 


47 (2) 


SGRQ total, units 


55.2 (2.0) 


CRQ total, units 


4.00 (0.13) 


VAS total, mm 


147 (9) 


Sputum total cell count X 10^^ cells/g. 


3.8 (2.8-5.3) 


geometric mean (95% CI) 




Sputum neutrophils, % 


69 (3) 


Sputum eosinophils, %, geometric 


1.3 (0.9-1.9) 


mean (95% CI) 





Data are presented as mean (SEM) unless otherwise indicated. 
CRQ = Chronic Respiratory Disease Interviewer- Administered Ques- 
tionnaire; ICS = inhaled corticosteroid; SGRQ = St. George Respiratory 
Questionnaire; VAS = visual analog scale. 
^Postbronchodilator. 



identify the presence of EM was 0.87 (95% CI, 
0.79-0.96) (Fig 2). All of the conventional measures 
of lung function (FEV^, DlcoA^a, RV, and FEVi/FVC) 
were reliable in distinguishing between EM only 
and EM plus BE and/or BWT (EM + BE/BWT), 
with receiver operator characteristic curves of < 0.6 
(Fig 3). 

In contrast to this, patients with radiologic EM in 
the presence or absence of BE and/or BWT did not 
have any significant difference in markers of proximal 
airway inflammation (sputum total cell count and cel- 
lular differential count) or bacterial load (measured 
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Figure 1. Nonproportional Venn diagram representing patients 
with EM, BE, and BWT on CT scanning. BE = bronchiectasis; 
BWT = bronchial wall thickening; EM = emphysema. 
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Table 2 — Clinical Comparisons in Patients With Demonstrable Absence or Presence of CT Scan Evidence 

of EM, BE, or BWT 



Characteristic 


Absence of EM/BE/BWT (n = 13) 


Presence of EM/BE/BWT (n = 62) 


P Value 


Male, No. (%) 


8(62) 


50 (81) 


0.63 


Age, y, mean (range) 


68 (49-88) 


67 (43-85) 


.90 


Pack-y smoked mean (range) 


35 (10-68) 


50 (10-153) 


.10 


Exacerbation history in previous y 


3.7 (0.9) 


4.4 (0.3) 


.43 


ICS dosage, |jLg 


831 (213) 


1548 (95) 


<.001 


FEVi % predicted^ 


62 (5) 


43 (2) 


<.001 


FEV/FVC,=^ % 


59 (3) 


45 (2) 


<.001 


SGRQ total, units 


50.3 (5.9) 


56.4 (2.0) 


.23 


CRQ total, units 


4.21 (0.35) 


3.95 (0.14) 


.45 


VAS total, mm 


149 (16) 


146 (10) 


.92 


Sputum total cell count X 10^ cells/g. 


3.5(1.9 -6.6) 


3.9 (2.7-5.6) 


.82 


geometric mean (95% CI) 








Sputum neutrophils, % 


66 (5) 


70 (3) 


.64 


Sputum eosinophils, %, 


1.2 (0.5-2.8) 


1.3 (0.9-2.0) 


.83 



geometric mean (95% CI) 



Data are presented as mean (SEM) unless otherwise indicated. BE = bronchiectasis; BWT = bronchial wall thickening; EM = emphysema. 

See Table 1 legend for expansion of other abbreviations. 

'iPostbronchodilator. 



by colony-forming units and abundance of 16S ribo- 
somal subunit encoding genes) (Table 3). 

Static Lung Volume and Transfer Factor 
Are Independent Predictors of EM in 
COPD, But Not of BE or BWT 

The DlcoA^a and RV were confirmed by logistic 
regression analysis to be independent predictors of 
EM on CT scanning in patients with COPD. In con- 
trast, the variables entered into the model were not 
independent predictors of BE or BWT (Table 4). 




Figure 2. Receiver operator characteristic curve (sensitivity and 
1-specificity) of diffusing capacity of the lung for carbon monox- 
ide corrected for alveolar volume for identifying emphysema. 



Conventional Lung Function Fails to Identify 
Patients With BE and BWT in COPD 

Having confirmed that static lung volumes and 
transfer factor were predictive of radiologic EM, we 
sought to identify whether these markers could inde- 
pendently group patients with EM -predominant 
disease (EM + ) or disease without EM (EM — ) in the 
presence or absence of BWT/BE, using unsupervised 
statistical techniques. Principal component analysis 
using airway physiologic parameters in all patients with 
full physiologic testing and CT scan (n = 64) identi- 
fied three factors with eigen values > 1 (Table 5), 
highlighting the following dominant physiologic com- 
ponents: (1) RV, (2) DlcoA^a, and (3) lung capacity 
Hierarchic cluster analysis determined three cluster 
groups. Two were EM predominant, discrete only 
in their degree of air trapping and gas transfer. The 
third cluster was characterized by a heterogeneous 
combination of EM and BE, with preservation of gas 
transfer and lung volumes (Table 6). 

Repeatability of Airway Physiologic Phenotypes 

Forty-seven patients had spirometry performed 
12 months after the initial characterization visit. In 
these patients, the repeatability of the measures of 
FEV^ and FVC was good, with intraclass coefficients of 
0.91 and 0.82, respectively (F< .001). There was no dif- 
ference in the mean change in lung function for each 
of the physiologic clusters identified after 12 months 
(cluster 1 [EM-predominant, moderate lung func- 
tion], 0.03; 95% CI, -0.11 to 0.17; cluster 2 [heteroge- 
neous, lung function preserved], —0.13; 95% CI, —0.23 
to —0.03; and cluster 3 [EM-predominant, severe lung 
function], 0.00; 95% CI, -0.12 to 0.12; P = .17). 
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with COPD, whereas the presence of BWT alone is 
sparse. Current guidelines use bedside spirometry and 
FEV^ to diagnose and guide the severity of COPD.^ 
Although this tool can be widely applied, it has less 
accuracy when used to measure small airways dys- 
function (airways with internal diameter < 2 mm).^ 
Using CT imaging, we have shown that, in patients 
with COPD, there are variation and overlap in the 
associated pathologic findings. Thus, the use of FEV^ 
alone to diagnose and guide management in COPD 
highlights that spirometry cannot define these patho- 
logic processes with clarity. Our finding of BE in those 
with spirometric classification of COPD is compara- 
ble to those observed previously^^ and highlights 
the heterogeneity that exists, which can, in turn, be 
partly delineated by radiologic characterization. ^s-^o 

In patients with EM on CT scan, we found wors- 
ened lung function, airflow obstruction, and larger 
increases in residual lung volume, compared with 
patients with radiologic evidence of BE or BWT, 
although the clinical parameters of airway inflam- 
mation, health status, and microbiologic qualitative 
or quantitative inflammation are indistinguishable 
between these radiologic phenotypes. This is in keep- 
ing with previous physiologic observations^^ and 



Table 3 — Airway Inflammometry^ Physiology^ and Health Status in Patients With EM Alone, 
BE and/or BWT Without EM, and EM With BWT and/or BE 



Characteristic 


EM (n = 29) 


BE/BWT-EM(n=12) 


EM + BWT/BE (n = 21) 


? Value 


Male, No. (%) 


22 (76) 


9(75) 


19 (91) 


.37 


Age, y, mean (range) 


64 (43-83) 


71 (60-84) 


69 (47-85) 


.07 


Ex-smoker, No. (%) 


17 (59) 


11 (92) 


12 (57) 


.09 


Pack-y smoked, mean (range) 


48(18-120) 


47 (10-153) 


55 (10-134) 


.68 


Exacerbation history in previous y 


4.0 (0.5) 


3.9 (0.9) 


5.1 (0.5) 


.32 


BMI, kg/m2 


25.2 (0.9) 


28.2(1.4) 


24.7 (0.4) 


.10 


Sputum total cell count X 10^ cells/g. 


3.2 (2.0-5.2) 


8.4 (4.3-16.5) 


3.6 (1.7-7.7) 


.18 


geometric mean (95% CI) 










Sputum neutrophils, % 


70 (4) 


77 (6) 


65 (6) 


.44 


Sputum eosinophil, %, geometric mean (95% CI) 


1.6 (0.9-2.8) 


0.8 (0.4-1.7) 


1.3 (0.6-2.7) 


.55 


Potential pathogenic organism, % 


14 


25 


38 


.08 


Sputum neutrophil count X 10^ cells/g. 


2.3(1.3-4.0) 


6.3 (2.8-14.1) 


1.9 (0.8-4.8) 


.18 


geometric mean (95% CI) 










Total bacterial load (16S), geometric mean (95% CI) 


1.58 (6.6"-3.28) 


1.78(3.5"-8.P) 


1.58 (3.2"-6.78) 


.99 


Colony forming units per mL of sputum, 


1.46 (6.35-3.06) 


1.26(1.56-9.16) 


7.85 (3.05-2.06) 


.71 


geometric mean (95% CI) 










FEV^ % predicted^ 


41 (3) 


59 (5) 


37 (4) 


<.01 


FEV/FVC,- % 


42 (2) 


56 (4) 


44 (2) 


.01 


Reversibility, % 


4.0(1.9) 


3.0 (2.9) 


-0.4 (3.3) 


.44 


RV, % 


162(10) 


115 (10) 


139 (10) 


.02 


Total lung capacity, % 


124(10) 


99 (5) 


103 (7) 


.13 


Dlco, mL CO/min/mm Hg 


4.0 (0.3) 


5.8 (0.6) 


3.6 (0.3) 


<.01 


Dlco % predicted 


46 (4) 


73 (7) 


32 (4) 


<.01 


DlcoA^a, mL CO/min/mm Hg/L 


0.8 (0.1) 


1.0 (0.1) 


1.3 (0.4) 


.37 


Dlco/Va % predicted 


56 (5) 


89 (10) 


49 (4) 


<.01 


SGRQ total, units 


57.9 (3.1) 


52.3 (3.5) 


56.4 (3.6) 


.61 


CRQ total, units 


3.79 (0.21) 


4.07 (0.31) 


4.11 (0.26) 


.57 


VAS total, mm 


163 (13) 


143 (25) 


126 (18) 


.24 



Data are presented as mean (SEM) unless otherwise indicated. Dlco = diffusing capacity of the lung for carbon monoxide; Va = alveolar volume; 

RV = residual volume. See Table 1 and 2 legends for expansion of other abbreviations. 

aPostbronchodilator. 




00 0^ 04 06 OS 10 
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Figure 3. Receiver operator characteristic curves for distin- 
guishing the presence of emphysema alone vs emphysema plus 
bronchiectasis and/or bronchial wall thickening. Dlco/Va = 
diffusing capacity of the lung for carbon monoxide/alveolar volume; 
RV = residual volume. 

Discussion 

In this study, we have shown that there is overlap 
in radiologic evidence of EM, BE, and BWT in patients 
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Table 5 — Extracted Factors With Eigenvalues > 1 Using 
Unsupervised Principal Component Analysis for Airway 
Physiologic Parameters 



Correlation Matrix 




Factors 




1 

1 


2 


1 

3 


FEV/FVC, % 


-0.49 


0.60 




FEV^ % predicted 


-0.35 


0.44 


0.74 


VC, % 






0.94 


FRC, % 


0.89 


-0.40 




RV, % 


0.85 






TLC, % 


0.90 






Dlco, % 




0.93 




DlcoA^a, % 




0.95 




FRC = functional residual capacity; TLC = total lung capacity; 
VC = vital capacity. See Table 3 for expansion of other abbreviations. 



evidence of a lack of correlation between sputum 
neutrophils and EM score in patients with COPD/^^i 
and outlines the radiologic complement to pheno- 
typing COPD. The evidence of eosinophilic airway 
inflammation and CT imaging in patients with COPD 
has been conflicting.^^'^^ However, these differences 
could be accounted for by differences in study design. 
Interestingly, Miller and colleagues^^ found close 
correlations between eosinophil degradation prod- 
ucts and EM score in a group of well-characterized 
patients with COPD, suggesting that these may con- 
tribute to the progression of EM. 

In this study, we have shown that FEV^ using 
regression model analysis was not predictive for the 
identification of EM, BE, or BWT on CT scan. This 
further illustrates the reduced accuracy of FEV^ in 
identifying the key pathologic processes of COPD. 
^ On the other hand, CT scan COPD phenotype 

identification has been associated with key clinical 
outcomes, including bronchodilator reversibility,^^^^ 
effect of airflow obstruction,^! and severity of COPD 
^ exacerbations. 44 We also concluded that the use of 

CD 

(^ conventional lung function testing alone was unable 

^ to delineate between patients with EM alone and 

g, those with EM plus phenotypes. 

Using cluster analysis, an unbiased statistical 
approach investigating airway physiologic pattern 

M expression, we identified three clusters that were 

^ ^ largely clinically indistinguishable in the expression 

S ^ of airway inflammation, bacterial load, and quality- 

^ of-life scores. These clusters identified two emphyse- 

I matous groups discrete by the degree and severity of 

^ hyperinflation and gas exchange, and a third group that 

g ^ § was heterogeneous. Although largely useful, physiology 

^ ^ alone could not clearly and robustly identify the patho- 

^ " ^ I So^ logic processes seen in COPD. The application of 

^ ^ S ^1 I CT imaging provided additional clinically important 

^ ^ ^ ^ ^ g information as a noninvasive biomarker in COPD.^o 

> 8 ^ 3 X I 1^ One limitation of this study is that we used qualita- 

Q pc; e2 pq tive and not quantitative classification of EM, BE, and 
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Table 6 — Demographic and Clinical Characteristics of Physiologic Clusters in COPD 



Cluster 1 (EM-Predominant, Cluster 2 (Heterogeneous, Cluster 3 (EM-Predominant, 
Moderate Lung Function Lung Function Preserved Severe Lung Function 
Characteristic [n = 30]) [n = 20]) [n = 14]) F Value 



Male, i\io. [ /c ) 


9Q /77\ 
Zo \ l 1 1 


1 7 {^K\ 
i / \oD) 


1 9 l9.f\\ 
IZ (,oD/ 


.Do 


Age, y, mean (range) 


Do 


70 


PA 
D4 


9Q 
.Zo 


Pack-y smoked, mean (range) 


Do 


Q9 
oZ 


AA 
44 


^ Ol 
<- .Ul 


Exacerbation history in previous y 


4.0 \\J.O) 


o.o iU.Dj 


A P (C\ P\ 

4.D [yj.'o) 


9K 
.ZD 


jjivii, Kg/m 


Z4.D \y).o) 


Zo.O il.Zj 


9 A Q /T o\ 

Z4.y ii.oj 


m 


Sputum total cell count X 10^ cells/g. 


4.Z (^o.o-D.oj 


4. / \o.D-o.Uj 


0.4 \A.l-o.V)) 


7Q 


geometric mean (95% CI) 










Sputum neutrophils, % 


P.^ (K\ 
DO \0) 


Do i_0/ 


7Q /fi^ 

/J (d; 


77 


Sputum eosinophils, %, 


l.D \L.Zj-Zj.\) ) 


U.O l^U.D-i.4/ 


9 0/1 Q zt 
Z.U i^l.0-4.0/ 


1 7 


geometric mean (95% CI) 










Potential pathogenic organism, % (95% CI) 


AA ('9Q fil ^ 


11 (^Z-o4j 


90 IK K9\ 
ZU \0-OZj 


HQ 

.uy 


Absolute neutrophil count X 10^ cells/g. 


Z.o {Z.yj-D.yj) 


Q fk /9 9 i=C C\ 
O.U [Z.Z-D.O) 


9 Q /I A P P\ 
Z.o (_l.4-D.Dj 


.yl 


geometric mean (95% CI) 










Total bacterial load (16S), 


Z.D l^l.D -D.D^ 


1 48 /7 Q7 A 98\ 
1.4 / .y -4.Z / 


D.O l^O.D -Z.4"/ 


.0\J 


geometric mean (95% CI) 










Colony- forming units. 


7 95 (A ^5 1 QS") 


1 ^6 (Q A5 Q P,6) 
l.D l^y.4 -O.D / 


9 Q6 M 76 Q 96\ 

z.y 1. 1 -o.Zi / 


1 7 
.1 < 


geometric mean (95% CI) 










FFV r>rf='rlir'f<=>r|a 
r j-j V 2^ /o uicU-icLeu. 




62 (3) 


97 ("^) 

Zj I \0) 


< 01 


r III Vj^/r Vvj," /c 


A^ (9) 

40 (,z; 


Do \0) 


'XK (9\ 
OD 


^ 01 
\ .Ul 


Jievei aiuiiiLy, /o 




D {OJ 




.12 


JAV, /C 




OO \\j I 


1 Q4 / 1 1 


^ 01 
\ .Ul 


TT n CK, 


1 1 n /4^ 


OD \D) 


1 94 (A\ 
1Z4 1^4/ 


^ 01 
\ .Ul 


lyLiCU/ V A /o preuicieQ 


fi*^ (A) 

DO (,4; 




49 (A^ 
4Z (^4; 


•cT 01 
\ .Ul 


o^^rvv^ Lorai, units 


04 [^) 


DO \Dj 


Di (,4; 


7Q 

. /y 


CRQ total, units 


3.8 (0.2) 


4.3 (0.3) 


4.3 (0.3) 


.30 


VAS total, units 


154 (13) 


164 (16) 


102 (18) 


.03 


VAS dyspnea, mm 


53 (4) 


44 (6) 


38 (8) 


.13 


VAS cough, mm 


41 (5) 


45 (6) 


25 (8) 


.12 


VAS sputum production, mm 


35 (4) 


38 (7) 


28 (8) 


.53 


VAS sputum purulence, mm 


25 (4) 


39 (6) 


11(3) 


<.01 


EM, % (95% CI) 


80 (62-91) 


25 (11-47) 


93 (66-100) 


<.01 


BE, % (95% CI) 


23 (12-41) 


35 (18-57) 


21 (7-48) 


.58 


BWT, % (95% CI) 


37 (22-55) 


5 (0-25) 


29 (11-55) 


.04 



Data are presented as mean (SEM) unless otherwise indicated. See Table 1, 2, 3, and 5 legends for expansion of abbreviations. 
'^Postbronchodilator. 



BWT. Although this qualitative description of EM, BE, 
and BWT introduces bias into a study such as this 
one, we aimed to reduce this by showing substantial 
agreement between two observers in using this form 
of radiologic description. Future studies that use 
similar unbiased statistical approaches combining mul- 
tiple dimensions of airway inflammation, physiology, 
and structure determined by quantitative densitom- 
etry and geometry are required. The study investigated 
patients who were concurrently enrolled in a longitu- 
dinal "biomarkers in COPD exacerbation" study. This 
may have introduced selection bias in the physiologic 
analysis because the indication for CT scan was pri- 
marily for clinical purposes, namely malignancy, BE, 
and lung volume reduction surgery. However, as com- 
monly found in clinical practice, the group with COPD 
was likely to have measures of lung function with 
imaging in a selected group only; our study thus 
highlights that lung function tools alone do not 
clearly identify important pathologic processes that 
may subsequently alter management. The time inter- 



val between characterization and CT scanning had 
a median of 15 months and may have affected the 
stability of the physiologic phenotype. We found in a 
large cohort of subjects with COPD that parameters of 
spirometry (namely, FEV^ and FVC) were unchanged 
12 months after characterization, and that measures 
of detailed lung function were also repeatable in this 
cohort. This further confirms our finding that the 
use of physiology that is repeatable can identify an 
EM phenotype, whereas imaging is required to iden- 
tify BE and BWT in COPD. Whether the presence of 
radiologic BE and BWT in COPD alters the clini- 
cal prognosis, management, and therapeutic response 
needs to be examined further in larger clinical stud- 
ies with specific phenotypic strategies of trial design. 

Conclusions 

In conclusion, COPD is a heterogeneous disease 
currently diagnosed and classified according to spirom- 
etry. The application of radiologic imaging can be a 
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useful tool in further "phenotyping" patients with 
COPD where clinical markers of airway inflamma- 
tion and health status are indistinguishable, and where 
static or dynamic lung function tests provide some, 
but not all, clarification. 
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